A  COMPLIANCE  K  CALIBRATICti  FOR  A  PRESSURIZED 


THICK-WALL  CYLINDER  WITH  A  RADIAL  CRACK 
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The  K  calibration  for  an  internally  pressurised, 
thick-wall  cylinder  with  a  straight,  radial  notch  h*« 
been  determined  from  a  compliance  test.  The  Bathed 
suggested  by  Irvin  is  used  wit'  compliance  defined  as 
tl4«  change  in  internal  volume  of  a  cylinder  divided 
by  applied  hydrostatic  pressure  rather  than  the 
usual  load-elongation  definition.  The  derivative  of 
internal  volume  change  with  respect  to  notch  depth, 

”a%  is  obtained  by  numerical  analysis  of  tangential 
strain  measurements  on  the  OD  of  the  test  cylinder. 

Thia  derivative  leads  directly  to  the  K  calibration 
for  the  cylinder.  Cubic  spline  functions  are  used 
to  approximate  both  the  strain  as  a  function  of 
posit in  on  the  cylinder  and  the  resulting  volume 
change  as  a  function  of  "a" .  Also  included  in  the 
determination  of  K  is  a  proof,  using  the  divergence 
theorem  in  the  theory  of  elasticity,  that  the 
derivatives  with  respect  to  "a"  of  internal  and 
external  volume  change  are  identical.  This  allows 
the  use  of  external  strain  measurements  to  determine 
K  bp «ed  on  internal  volume  change. 

The  compliance  K  calibration  nearly  coincides  with 
a  semi-infinite  plate  solution  simulating  both  the 
tangential  stress  due  to  pressure  and  the  direct 
effect  of  pressure  in  the  notch,  Kj  ■=  1.12  <f$  *f~W a  + 
1,13  p  //rr  al  Thia  unexpected  agreement,  particularly 
for  values  of  a/w  up  to  0.6,  is  explained  by  the 
combination  of  bending  constraint  and  drop  off  of 
tangential  stress  in  ths  cylinder  wall. 
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Crack  or  notch  depth 

Cro>,3-8ectlon  area  of  cylinder  material 

Area  enclosed  by  the  outside  perimeter  of  the  cylinder 

Compliance,  change  in  length  divided  by  load  for  a  unit  thickness 
sheet 

Compliance,  internal  volume  change  divided  by  pressure  for  a  un 
length  cylinder 

Young's  modulus 

Crack  extension  force  per  unit  length  along  the  crack  front 

Constant  applied  to  plastic  zone  correction 

Opening  mode  stress  intensity  factor 

Hydrostatic  pressure  on  the  ID  of  a  cylinder 

Radial  coordinate  of  a  cylinder 

Inner  radius  of  a  cylinder 

Outer  radius  of  a  cylinder 

Volume  enclosed  by  tin  inner  surface  of  a  unit  length  cylinder 
Material  volume  of  a  unit  lengtn  cylinder 

Volume  enclosed  by  the  outer  surface  of  a  unit  length  cylinder 

Wall  thickness  of  a  cylinder 

Length  coordinate  of  „  cylinder 

Uniform  tensile  stress 

Tangential  stress  in  a  cylinder  wall 

Tensile  yield  stress,  0.25t 


Other  notation  defined  in  the  text. 
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INTRODUCTION 


The  determination  of  stress  intensity  factors  by  an  experimental  compli¬ 
ance  test  is  seldom  reported  in  the  literature.  This  is  probably  due  to  the 
difficulty  in  performing  compliance  experiment*  with  suffic*  -nt  accuracy  and 
also  due  to  the  limited  application  of  the  experimental  results,  once  obtained. 
Numerical  and  analytical  X  calibrations  are  often  at  least  as  reliable  and  as 
easy  to  obtain.  They  can  often  be  applied  to  many  different  geometries  with 
little  difficulty.  However,  experiment  ad  K  calibrations  cam  be  useful  when 
the  geometry  is  difficult  to  model  or  when  the  experiment  is  relatively 
simple. 

In  this  case  the  compliance  testing  of  a  pressurized  thick-wall  cylinder 
with  a  narrow  radial  notch  was  performed  as  part  of  another  experiment.  Our 
purpose  here  is  to  use  the  available  compliance  results  to  determine  the  K 
calibration  of  a  thick-wall  cylinder  for  use  in  design  and  J.or  comnarison 
with  other  fracture  mechanics  analyses.  The  analysis  in  the  literature  which 
appears  close  t  to  a  thick-wall  cylinder  geometry  is  Bowie's  solution  for  a 
pressurized  hole  with  a  radial  crack  in  an  infinite  plate^.  His  solution 
should  correspond  to  a  pressurized  cylinder  with  a  radial  crack  in  an 
infinitely  thick  wall.  The  experimental  results  will  be  compared  with  tills 
and  other  analyses  in  the  sections  that  follow. 

FRACTURE  MECHANICS  ANALYSIS 

The  compliance  method  used  to  determine  the  X  calibration  for  a  pressur¬ 
ized  cylinder  is  patterned  after  Irwin's  method^ .  The  usual  definition  of 
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compliance,  l.e.,  the  ratio  of  change  in  length  of  a  specimen  to  the  applied 
1  ad,  is  replaced  by  the  ratio  of  change  in  internal  volume  of  a  cylinder  to 
the  applied  pressure.  Internal  volume  and  pressure  are  used  to  obtain  an 
expression  for  crack  extension  force  for  a  pressurized  cylinder  similar  to 
Irwin's  expression  for  a  load-displacement  system  from  Ref.  (2),  which  is: 

G  =  7  r2£ 

a  i 

where  compliance,  C  =  -jp  ,  is  in  terms  of  the  tensile  load,  F,  on  a  specimen 
and  the  resulting  change  in  length,  Al.  In  the  above,  G  and  C  are  written 
in  terms  of  unit  specimen  thickness,  i.e.,  unit  dimension  along  the  crack 
front. 

For  a  unit  length  cylinder  at  constant  internal  pressure,  the  differen¬ 
tial  change  in  strain  energy,  U,  corresponding  to  a  change  in  internal  volume 
can  be  written 


dU  =  |  p  d(  AVt) 

The  Griffith  crack  extension  force  is  defined  as  the  differential  change  in 
strain  energy  of  a  system  due  to  a  change  in  crack  depth.  Thus, 

dU  =  r  =  1  r  dCAVi) 
da  2  3a 


Uwing  a  modified  definition  of  compliance  Cv 
the  crack  extension  force  becomes 


and  with  p  constant, 


-  p2 

2  -a; 


(2) 
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an  expression  analogous  to  eq.  (1)  with  tensile  load  replaced  by  hydrostatic 
pressure.  G  and  Cy  are  written  in  terms  of  unit  dimension  along  the  crack 
front  as  before,  which  in  this  case  is  in  the  length  or  z  direction  of  the 
cylinder,  see  Fig.  1.  For  plane  strain  conditions,  as  in  the  compliance  test 
described  below,  the  stress  intensity  factor  is 


1  -  P 


S  P2  d(AV±/p) 
2(1  -  P  2)  da 


The  K  calibration  can  be  calculated  simply  by  determining  the  change  in 
internal  volume  cf  the  cylinder  at  a  given  value  of  pressure  for  a  series  of 
different  notch  depths  and  then  performing  the  indicated  differentiation. 
However,  the  measurement  which  can  be  easily  made  is  the  change  in  outside 
volume  of  the  cylinder.  If  it  can  be  shown  that  the  derivative  of  compliance 
with  respect  to  crack  depth  determined  from  the  outside  volume  change  is 
equivalent  to  that  determined  from  the  inside  volume  change,  then  the  problem 
is  solved.  Rice^)  has  outlined  an  analysis  which  indicates  that  the  deriva¬ 
tives  are  equivalent.  In  the  complete  analysis  below  we  prove  that  the 
change  in  material  volume  of  a  cylinder  is  unaffected  by  the  presence  of  the 
notch  provided  that  tractions  on  the  opposite  sides  of  the  notch  are  the  same, 
and  thus  the  two  derivatives  of  compliance  are  equivalent. 

The  change  in  material  volume  in  an  elastic  body  is  given  by 


'»  ■/ 


For  plane  strain,  the  dilatation  e  -  £  i  "1,2,  hence 


A  V*/unit  length 
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r:s*K,  i  ii  WB w* 


where  A  la  the  cross-aecticn  area  of  tfte  body. 

The  stress-strain  relations  for  plane-strain  are 


contraction  gives 


Substituting  (5)  into  (4)  gives 


A  V^/uni'c  length 


(5) 


(6) 


The  equation  of  equilibrium,  in  absence  of  body  force#  is 


Tik,k  "  0  X'k  =1'2 

Now  consider  the  following  integral  with  the  help  of  the  above 

/  TiJ  ^  *  f  Tik  ^kj  ^  +  /  Xj  Tik,k 

A  A  A 

=  /  (  rik  xj,k  +  Tik,k^  ^ 

A 

“  /  (  Tik  XjJ,k 

A 

Using  the  divergence  theorem  we  get  the  following  line  integral 

I  (  rik  x^,k  &  *  J  Tik  XJ  ^  “  /  Ji  Xj  ^ 

A  S  S 
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where  are  the  ^^ectiori  cosines  of  the  exterior  normals  £  to  the  line 

»  )> 

S,  ar.d  Tj  are  the  components  of  the  tractions  JT  acting  on  S,  see  Ref.  (4), 

hence  we  have 


/  T«  <*  -  \  h 


XJ  ** 


an*..  contraction  gives 


J  T. .  <*  .  f  Xj  dS  -  /  T  •  r 


Substituting  (?)  into  (6)  we  have 


4Vm/unit  length  = 


T  •  r  dS 


Breaking  the  line  into  four  parts,  S  =  +  S2  +  +  S^s  which  correspond 
to  the  cylinder  and  notch  surfaces  as  shown  in  Fig.  1,  and  keeping  the  outer 


normals  in  mind,  we  get 


AVjj/unit  length  = 


0  Tx  .  r  -  |)  T2  •  r  > 

>1  S2  \ 

Fv^a  )  ; 


X3  •  r  dS 


T,  •  r  dS 
4 


v  V  £  tr 

where  T-, ,  T~  are  tractions  on  the  inner  and  outer  radii  and  T^,  17  are  th« 

tractions  on  the  bottom  and  top  of  the  notch.  Now  when  To  *  $l,  i.e.,  when 

the  tractions  on  both  sides  of  the  notch  are  the  same,  the  last  two  integrals 

vanish,  =nd  the  proposition  is  proved.  This  is  the  case  in  the  problem  under 

consideration. 
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For  our  cylinder  then,  we  have  f,  »  „  ®_  p  To  *  0,  where  e_  is  a  uru.i 

•\*X  wi  * 

vector  in  the  radial  direction.  Using  eq.  (8)  we  get 


dVy/u nit  length 


the  change  in  material  volume  of  the  cylinder  per  unit  length  in  the  s 
direction.  This  result  can  be  checked  directly  from  the  Lame  solution  (5). 

From  the  above  we  have  yj  ( A  Vm/unlt  length)  *0  so  that  eq.  (3) 
can  be  written  in  terms  of  the  outside  volume,  V0,  of  the  cylinder 

E  d{  A  V^/p) 

2(1  -  da 

COMPLIANCE  TEST 

A  1.028  in.  ID,  1.995  in.  CD,  6  in.  long,  4340  ateel  cylinder  was 
pressure  tested  as  shown  in  Fig.  2.  The  cylinder  was  tested  in  an  open  ended 
condition,  i.e.,  the  rams  shown  in  the  sketch  served  carry  the  end  loads 
due  to  the  pressure  as  well  as  seal  the  ends  of  the  cylinder.  Thus,  there 
were  no  *  direction  stresses  produced  in  the  cylinder  due  to  pressure  end¬ 
loads.  The  pressure  fluid  used  was  a  synthetic  instrument  oil.  No  attempt 
was  made  to  seal  the  notch  area,  so  effects  due  to  pressure  in  the  crack  were 
present  in  the  tube. 

The  pressure  was  applied  in  four  10,000  psi  increments  for  each  of  ten 
notch  depths  from  0  to  .316  inches.  The  .015  in.  thick,  4  in.  long  notch 
was  cut  deeper  following  each  pressure  run  by  an  electrical  discharge  machin¬ 
ing  process.  The  circumferential  strain  on  the  CD  of  the  cylinder  was 
measured  at  the  center  of  notch  length  with  16  resistance  strain  gages  around 
the  tube.  The  strain  gage  data  for  10,000  psi  pressure  are  listed  in  Table  I. 
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Figure  2.  Caapliance  Test  Set-up 


Table  I.  Measured  Tangential  Strain  at  10,000  pel  Pressure 


'’AGEE  LOCATION  _ _  NOTCH  LENGTH,  In. 


WO. 

Degrees 

0 

*043- 

£$1 

.II? 

.152 

■1  — Sinii  1 

.181 

.189 

,330 

.3lS 

1 

0 

240 

220 

145 

80 

5 

-  10 

-  90 

-120 

-375 

-495 

2 

+  7 

230 

230 

ISO 

145 

85 

95 

90 

235 

450 

3 

-  7 

240 

225 

160 

100 

105 

30 

-  30 

-  5 

-  40 

100 

4 

+  14 

240 

255 

235 

265 

265 

285 

320 

325 

675 

135 

5 

-  14 

235 

240 

215 

210 

230 

210 

225 

255 

555 

805 

6 

+  29 

240 

265 

310 

360 

405 

420 

470 

495 

760 

840 

? 

-  29 

240 

265 

300 

340 

400 

400 

455 

485 

780 

875 

8 

+  57 

240 

260 

285 

325 

365 

360 

400 

425 

580 

635 

9 

-  57 

235 

265 

290 

330 

365 

380 

405 

440 

620 

665 

10 

+  86 

220 

250 

260 

280 

295 

285 

295 

320 

395 

405 

11 

-  86 

230 

255 

225 

265 

— 

— 

315 

. — . 

— 

— 

12 

+115 

240 

240 

245 

245 

240 

240 

215 

240 

245 

220 

13 

-115 

240 

250 

240 

255 

250 

250 

240 

255 

— 

250 

14 

+143 

240 

240 

230 

215 

205 

200 

180 

190 

140 

no 

15 

-143 

240 

240 

— 

220 

205 

215 

190 

205 

155 

120 

16 

±180 

240 

240 

220 

210 

190 

185 

165 

160 

95 

45 

The  gage  positions  shown  are  in  degrees  away  from  the  point  on  the  CD  directly 
above  the  notch.  A  plot  of  the  variation  of  strain  around  the  tube  for  a 
deep  notch  and  30,000  pei  pressure  is  »hown  in  Fig.  3  compared  with  the 
uniform  strain  value  for  an  uncracked  cylinde.  at  the  same  pressure^).  The 
higher  average  strain  in  the  tube  due  to  the  presence  of  the  notch  can  be  seen 
by  a  simple  visual  comp...  tson  of  the  two  plots. 

By  careful  numerical  analysis  of  all  the  strain  gage  data,  the  change  in 
outside  volume  of  the  cylinder,  A  V0,  and  thus  its  compliance  can  ba  accurately 
determined  as  a  rune t ion  of  notch  length.  Since  the  strain  data  is  taken  at 
the  center  of  the  6  in.  length  of  the  cylinder,  any  end  effects  will  have 
vanished  and  a  plane  strain  analysis  is  appropriate. 
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NUMERICAL  DATA  ANALYSIS 

By  considering  only  a  cross  section  of  the  cylinder  through  the  strain 
gages,  the  numerical  problem  may  be  simplified  and  stated  as  follows.  Given 
the  tangential  strain  «t  selected  points  on  the  circumference  of  a  circle, 
we  wish  to  determine  numerically  the  shape  of  the  perimeter  of  the  distorted 
figure  and  from  this  the  increment  in  area.  Lacking  knowledge  of  tne  radial 
displacement,  the  problem  is  somewhat  ambiguous.  However,  if  we  can  interpolate 
the  tangential  strain  measurements  with  an  angular  function  and  assume  that 
the  departure  of  the  distorted  perimeter  from  the  original  circle  is  small, 
then  two  simplified  approaches  suggest  themselves. 

First,  we  can  integrate  the  tangential  strain  to  obtain  the  change  in 
length  of  the  perimeter 


4P 


/2V 

r 


ee{e)  , 


and  assuming  the  distorted  figure  remains  circular,  the  change  in  area  enclosed 
by  the  perimeter  is 


A  K 


ZP  AP  +  (AP)2 
4  TT 


As  an  alternati/e  approach  we  can  assume  that  any  cu-splacement  of  the 
perimeter  is  solely  radial,  that  is 


<r2-,> 


and  then  from 
£  (**2> 


1_  (r2,9  )  Uf  (**2*  * ) 

r2 


£ 

'X 


B  3 
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we  have  ur  (r2,  0)  ~  ^^(r2»^)  *  r2 

f2ir  r  r2  +  ur'r2, 8  ) 

4A0  =  J  J  r  dr  d 0 

o  r2 

Since  these  two  approaches  give  reasonable  bounds  on  the  amount  of 
distortion  of  a  circular  cross  ae  tion,  we  can  assume  that  the  degree  to 
which  they  approximate  each  other  is  an  indication  of  the  number  of  signif¬ 
icant  figures  in  our  approximation  of  the  increment  in  area. 

Our  choice  of  an  angular  function  was  an  interpolating,  periodic,  cubic 
spline.  This  function  was  chosen  on  grounds  of  experience,  intuition  and 
personal  interest;  however,  it  does  belong  to  the  class  of  functions  which 
approximate  minimum  strain  energy  and  thus  is  a  likely  candidate  for  inter¬ 
polating  (or  approximating)  strain  readings. 

We  give  here  a  definition  of  a  cubic  spline.  Thc*e  who  are  interested 
in  the  details  of  the  construction  and  manipulation  of  such  functions  are 
•sferred  to  refs.  6->. 

Definition  of  a  Cubic  Spline 

Given  an  interval  oL  s  x  s  (5  > 

a  mesh  on  the  interval 

=  Xo«x i  ...  **  0  , 

and  an  associated  set  of  ordinates 

Y  5  70>  71*  •  •  •* 

then  a  cubic  spline  satisfies 

SA  (T;  x)  6  C2  on  [•«,/*] 

3$  (Y;  x»)  ™  y^  (J  “  0,1,  ...,  N) 
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and  is  coincided  with  a  cubic  on  each  subinterval 

x3-ls  x  SXJ  0-1.2,  ....  N) 

If  in  addition 

S 4(p,(«t+)  =  <!•-  0,1,2) 

the  spline  is  said  to  be  periodic  with  period  (  p  -  oC  ) . 

If  alternatively 

s4  (I;  Xj)  =  yj  +  <J  -  0,1, 

with  the  subject  to  some  minimizing  constraint,  it  is  said  to  be  an 
approximating  rather  than  an  interpolating  spline. 

As  described  in  the  previous  section,  the  data  was  taken  at  16  positions 
for  10  crack  depths  and  4  pressures.  Of  the  640  possible  readings,  16  were 
missing  because  the  tube  ruptured  before  the  last  scheduled  measurement,  35 
were  missing  because  of  gage  failure,  3  were  apparent  gage  failures,  and  one 
was  Judged  a  transcription  error.  The  last  39  were  replaced  by  bi-quadratic 
interpolation. 

To  indicate  the  smooth  nature  of  the  strain  data  and  interpolating 
functions  and  the  obvious  nature  of  the  transcription  error,  we  show  fitted 
spline  functions  for  four  pressures  at  a  crack  depth  of  0.117  inches,  see 
Fig.  4*  The  datum  at  +  86°  and  30,000  pei  looks  like  au  outlier.  It  was 
recorded  as  945  JJ*  in/ln  and  quite  likely  was  really  845  in/ in,  the  value 
on  the  dashec  curve  in  Fig.  4. 

It  ie  also  apparent  upon  close  examination  of  the  data  in  Fig.  4  that 
the  minimum  strain  is  not  at  the  indexed  zero.  Examination  of  the  sectioned 
tube  confirmed  that  the  gages  were,  indeed,  displaced  a  snail  positive  angle 
from  the  true  location  of  the  simulated  crack;  however,  this  had  no  effect 
on  the  computation. 
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POStT/ON  PEL4T/VE  TO  NOTCH,  DEOP  EES 


Evaluation  of  the  change  in  area,  4A0,  using  tho  two  methods  described 
were  found  to  agree  with  each  other  to  four  figures .  These  results  are 
presented  in  Table  II. 


Table  II.  Change  in  Area  Enclosed  by  the  Outside  Perimeter  of  the  Cylinder 

dAcjin.^ 


a,  in. 

at  10.000  psi 

at  20.000  psi 

at  30,000  psi 

0.0 

0.001480 

0.003009 

0.004581 

0.006157 

0.043 

0.001558 

0.003130 

0.004713 

0.006193 

0.083 

0.001549 

0.003126 

0.004770 

0.006371 

0.117 

0.001638 

0.003315 

0.005061 

0.006663 

0.140 

0.0U1689 

0.003452 

0.005211 

0.006°5r 

0.152 

0.001685 

0.003446 

0.005271 

0.00',J2> 

0.181 

0. 003.707 

0.003559 

0.005458 

0.007289 

0.189 

0.001817 

0.003682 

0.005602 

0.007436 

0.280 

0.002235 

0.004610 

0.007005 

0.009555 

0.316 

0.002416 

0.004939 

0.007643 

— 

It 

was  now  desired  to 

evaluate  eq.  (9)  repeated  below 

E 

d(  21  Vo/p)  * 

1/2 

Kj/p  = 

2(1  - 

V 2)  da 

(9) 

where  dV0  is  the  change  in  outside  volume  of  the  cylinder  per 

unit  length 

and  thus 

is  numerically  equivalent  to  A AQ,  , 

and  also 

*  , 

— 

E 

d(  flV0/p) 

1/2 

%  /  P  - 

2(1  - 

V  ^)  da* 

(10) 

2  n 


Kl 


is  an  effective  crac.i  length  suggested  by 


where  a*  a  + 

Irwin ^  to  correct  for  effects  of  plastic  deformation  at  the  crack  tip. 

We  also  performed  computations  varying  the  correction  by  c.  constant,  h, 

2 


h 

Trr 


(-) 

\  Oyl 


(11) 


to  examine  the  effects  of  modifying  the  correction  procedure. 
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To  approximate  4V0(a)  we  used  a  nor.  periodic,  approximating,  cubic 
spline  with  an  algorithm  which  minimized 

r Sri  ^  | 2 

J  ( s4”(a) }  da  +  I  s*(ai}  ■  **  i } 

0 

where  X  is  introduced  to  allow  us  to  strike  a  balance  between  the  afficu 

smoothing  desired  versus  our  leh  to  respect  the  integrity  of  the  data. 

As  a  guide  to  the  choice  of  ^  m  consider  the  type  of  Kj  expressio 

common  for  infinite  bodies  Kj  =  (constant)  ard  see  that 

d2  (Kl) 

da^ 

is  strictly  monotonic  increasing.  Thus,  we  want  X  to  be  small  enough  s 
that  this  is  also  true  of  our  approximating  spline. 

On  the  other  hand  we  must  of  necessity  evaluate  a*  and  K*/p  by  an 
iterative  process  and  X  must  not  be  so  small  that  this  process  is  unota! 
A  value  of  X  =  100  proved  to  bo  a  reasonable  cor  remise  leading  to  rapi 
and  stable  convergence  of  Kj/p. 

RESULTS  AND  DISCUSSION 

The  K  calibration  results  from  the  numerical  ar alysis  of  the  strain 
(eq.  9)  are  shown  in  Fig.  5.  An  ideal  experiment  e  id  analysis  would  be 
expected  to  give  results  wh'th  all  lie  on  the  same  curve.  We  believe  td 
the  main  factor  which  accounts  for  the  variation  c  the  Kr/p  curves  is  p 
deformation  at  the  notch  tip.  Other  factors,  sue}  as  the  finite  notch  w 
and  thermal  extension  of  the  cylinder  due  to  pres  ure  fluid  heating,  cou 
affect  the  results.  But  plastic  deformation  is  a  pec ted  to  oc;ur  most  r* 
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Figure  5.  Compliance  K  Calibration  Results 
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for  combined  high  pressure  and  deep  notches,  the  same  conditions  for  which 
the  variation  of  Kj/p  is  the  greatest.  Also,  the  strain  readings  taken  at 
zero  pressure  after  each  run  suggest  th&t  plastic  deformation  is  occurring. 
Ths  average  residual  strain  for  all  readings  through  .181  in.  notch  depth  is 
3  ti  in/in.,  8  ja.  in/in.  for  the  .189  reading,  and  70  ytt  in/in  for  the  .280 
in.  readings.  The  tube  fractured  (at  38,500  psi)  before  the  last  pressure 
increment  with  a  .316  in.  notch  present.  The  increase  in  residual  strain  and 
eventual  fracture  with  increasing  notch  depth  ia  a  clear  indication  that 
significant  plastic  deformation  was  present  in  the  cylinder  for  deep  notches 
at  high  pressure. 

The  K  results  including  the  plastic  zone  correction  of  eq.  (10)  are 
shown  as  dashed  curves  in  the  figure.  The  large  shift  of  the  high  pressure- 
deep  notch  results  has  little  significance  since  the  excessive  plastic 
deformation  probably  present  in  these  results  makes  a  fracture  mechanics 
analysis  inappropriate.  However,  it  is  interesting  to  see  in  Fig.  6  that  a 
smaller  correction  of  h  -  .3  in  eq„  (11)  causes  all  the  deep  notch  data  to 
fall  on  the  same  curve  within  0.5$.  This  tends  to  support  our  belief  that 
plastic  deformation  around  the  notch  is  the  major  uncertainty  in  the  experi¬ 
ment  and  analysis. 

The  corrected  10,000  psi  results  from  Fig.  5,  considered  to  be  the  best, 
are  repeated  in  Fig.  7  along  with  various  analyses.  Results  are  shown  fron 

/i  \ 

Bowie ‘s  analysis'  '  for  a  hole  with  a  radial  crack  in  an  infinite  plate  under 
biaxial  tension.  As  suggested  in  ref.  (1),  the  plane-strain,  plate  solution 
can  be  used  for  an  internally  pressurized  cylinder  of  infinite  wall  thickness 
with  a  radial  crack  if  tensile  load  is  replaced  by  pressure  in  the  solution. 
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The  resulting  infinite  cylinder  K  calibration  coincides  with  the  well-known 


semi-infinite  plate  solution 


(10) 


Kx  =  1.12  6  "\frir& 


(12) 


for  shallow  notches,  as  would  be  expected.  The  semi-infinite  plate  solution 
ia  plotted  in  terns  of  pressure  by  using  the  Lame  relation^)  between 
pressure  and  tangential  stress  in  a  cylinder. 


rfe  -  p 


(r2/r;  +  1 

"_\2 


(V1 


- 1 


(13) 


where  r  =  r-^  +  a 


The  maximum  value  of  6 $  (corresponding  to  the  ID)  was  used  for  plotting 
the  plate  solution,  so  th®  higher  Kj/p  values  for  deep  notches  are  expected. 

The  K  calibration  for  a  pressurized  cylinder  should  include  the  direct 
effect  of  pressure  in  the  notch  as  well  as  the  effect  of  the  stresses  in  th© 
cylinder  which  are  present  with  no  notch.  Bueckner^ 
for  a  pressurized  notch  in  a  half  space. 


'^11)  obtained  a  solution 


Kj  =  1.13  P  V  ^  a  , 


(14) 


where  the  constant  in  the  expression  can  be  thought  of  as  a  free-surface 
correction  factor  applied  to  the  solution  for  a  pressurized  notch  in  an 
infinite  plate The  superposition  of  the  direct  pressure  effect,  eq. 
(14),  and  the  tangential  stress  effect  using  the  maximum  value  of  Gq  , 
eqs.  (12)  and  (13),  results  in 


Kt  =  3.06  p  ~\frr  a 


(15) 
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The  good  agreement  between  the  semi-infinite  plate  solution  with 
combined  pressure  and  stress  and  the  experimental  results  is  surprising, 
particularly  for  deep  notches.  It  suggests  that,  if  the  experiment  is  to  be 
believed,  effects  mist  be  present  in  the  experiment  which  reduce  the  expected 
increase  in  Kj  for  deep  notches  in  a  finite  size  specimen.  The  two  dominant 
factors  which  lower  the  Kj  of  the  cylinder  are  the  drop  off  of  tangential 
stress  through  the  wall  described  by  eq.  (13)  and  the  doubly  connected  nature 
of  a  hollow  cylinder  which  tends  tc  prevent  ending  in  the  wall. 

Equation  (16)  is  an  attempt  to  account  for  these  factors  by  including 
the  decrease  in  as  a  function  of  r  and  by  including  &  bending-constraint 
expression  which  limits  the  increase  in  Kj  for  deep  notches.  Our  choice  for 
an  expression  to  represent  the  bending  constraint  is  the  ratio  of  the  Kj  for 

a  constrained-end,  SEN  plated)  to  the  Kj  for  a  semi-infinite  plate. 

finite  plate, 

-  constrained  ends,  - 

KI  ref, (13) _ 

_  a  end-infinite  * 

L  I  plate,  eq.  (12)  -* 

Of  the  two  sets  of  results  from  ref.  (13),  we  chose  those  for  a  plate  of 
2.6  length-to-width  ratio  as  the  better  representation  of  the  bending 
constraint  in  the  tes'.  cylinder.  Although  the  effects  of  Gq  and  bending 
constraint  both  change  by  as  much  as  50#  for  deep  cracks,  they  counteract 
one  another  so  that  the  neb  change  ie  hardly  noticeable.  An  expanded  portion 
of  the  experimental  results  and  values  from  eqs.  (15)  and  (16)  show  Just  a 
few  percent  spread,  see  Fig.  8.  So  apparently,  the  smaller  increase  in  Kj  for 
deep  notches  in  a  cylinder  with  decreasing  stress  through  the  wall  and  bending 
constraint  results  in  a  K  calibration  for  the  teet  cylinder  not  much  different 
from  that  of  a  semi -infinite  plate. 


KI  = 
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It  is  difficult  to  comment  on  th*  accuracy  of  the  experimental  results 
since  no  solution  is  available  for  a  direct  comparison.  Major  effects  dJe 
to  notch-tip  plasticity  are  probably  avoided  by  using  the  low  pressure  data. 

The  width  of  the  notch  is  the  other  factor  which  could  have  a  major  effect 
on  the  results.  The  good  agreement  between  the  experimental  K  calibration 
for  shallow  notches  and  the  infinite  plate  solution  indicates  that  the  notch 
width  had  no  major  effect  since  any  effect  should  be  at  least  as  great  for 
shallow  notches.  Finally,  an  indication  that  the  experimental  K  calibration 
is  not  grossly  in  error  is  the  Kj  value  calculated  from  the  experimental 
curve  of  Fig.  7  for  the  38,500  psi  fracture  pressure  and  .316  in.  notch  length. 
The  value,  120,000  psi  Jin.,  is  near  the  critical  value  of  Kj  expected  for 
the  cylinder  material. 

Finally,  it  should  be  emphasized  that  the  experimental  K  calibration 
shown  in  Fig.  7  strictly  applies  only  to  a  cylinder  of  the  dimensions  given. 
However,  we  believe  it  can  be  extended  to  any  cylinder  with  the  same 
ratio  by  using  the  square- root  size  factor  cannon  to  most  K  calibrations  of 
finite  specimens.  Thus,  a  Kj/p versus  a/w  plot  of  the  result*  could  be 
used  for  any  cylinder  of  the  same  diameter  ratio. 
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straight,  radial  notch  has  been  determined  from  a  compliance  test.  The  method  sug¬ 
gested  by  Irwin  is  used  with  compliance  defined  as  the  change  in  internal  volume  or*  a 
cylinder  divided  by  applied  hydrostatic  pressure  rather  than  the  usual  load-elongation 
definition.  The  derivative  of  internal  volume  change  with  respect  to  notch  depth, 

"a",  is  obtained  by  numerical  analysis  of  tangential  strain  measurements  on  the  0D  of 
the  test  cylinder.  This  derivative  leads  directly  to  the  K  calibration  for  the 
cylinder.  Cubic  spline  functions  are  used  to  approximate  both  the  strain  as  a  function 
of  position  on  the  cyli  der  and  the  resuming  volume  change  as  a  function  of  "a”.  Also 
included  in  the  determination  of  K  is  a  proof,  using  the  divergence  theorem  in  the 
theory  of  elasticity,  that  the  dt  -ivatives  with  respect  to  "a"  of  internal  and  external 
volume  change  are  identical.  This  al.’ows  the  use  of  external  strain  measurements  to 
determine  K  based  on  internal  volume  change. 

The  compliance  K  calibre* i<i  nearly  coincides  with  a  semi-infinite  plate 
solution  siawlating  both  the  tangential  stress  due  i^>  pressure  and  the  direct  effect 
of  pressure  in  the  notch,  Kj  »  1.12  v"  -■  1  ♦  1 . 13  p  v'  ~  1.  This  unexpected  sgree- 
ment,  particularly  for  values  of  a/w  up  to  0.6,  is  explained  by  the  combination  of 
bending  constraint  and  drop  off  of  tangential  stress  iip  the  cylinder  wall. 
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